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O . Abstract 
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We extend a recent three-loop calculation of nuclear matter in the systematic frame- 
work of chiral perturbation theory to finite temperatures T. The contributions from one- 
and two-pion exchange diagrams which cause nuclear binding and saturation at T = 
are included for T > in the density and temperature dependent free energy per particle, 
F(p,T). The so-called anomalous 2-7r-exchange contribution A(p,T) (with no counterpart 
in the ground state energy density at T = 0) is consistently included. The calculated 
I/"") | pressure isotherms display the familiar first-order liquid-gas phase transition of isospin 

symmetric nuclear matter with a critical point at T c = 25.5 MeV and p c = 0.09 fm" -3 . The 
■ too high value of the critical temperature originates from the strong momentum depen- 

dence of the underlying single-particle potential U (p, kfo) near the Fermi-surface. We also 
consider pure neutron matter at T > in the same framework and find fair agreement 
with sophisticated many-body calculations for neutron densities p n < 0.2 fm -3 . 
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Ctf ■ The analysis of data from low-energy heavy-ion collisions in the regime of nuclear fragmenta- 

tion has lead to the picture that heated nuclear matter undergoes a first-order phase transition 
from a liquid-like state to a vaporized gas state [1], 0. This liquid-gas phase transition of 
isospin symmetric nuclear matter is in fact very similar to that of the familiar van-der-Waals 
gas, with a corresponding critical temperature of T c ~ (16 — 18) MeV |2| and a critical density 
of p c = (0.06 — 0.07) fm -3 Clearly, the dynamical description of this phase transition is an 
important topic in any microscopic calculation of nuclear matter. In the <7u;-mean field model of 
Serot and Walecka || nucleons are described as Dirac-quasiparticles moving in self-consistently 
generated scalar and vector mean fields and a critical temperature of T c ~ 19 MeV is typically 
found [|]]. As an other example, the sophisticated many-body calculations of the Urbana group 
H using the V14 effective NN-interaction (plus an adjustable three-nucleon interaction) predict 
a critical temperature of T c ~ 18 MeV 0. For earlier work on this topic using effective Skyrme 
forces, see ref.Jj). 

In a recent work f?|], we have used chiral perturbation theory for a systematic treatment 
of the nuclear matter many-body problem. In this calculation the contributions to the energy 
per particle, E(kf), originate exclusively from one- and two-pion exchange between nucleons 
and they are ordered in powers of the Fermi momentum kf (modulo functions of kf/m n ). It 
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has been demonstrated in ref.H that the empirical saturation point and the nuclear matter 
compressibility K ~ 255 MeV can be well reproduced at order 0(kj) in the chiral expansion 
with just one single momentum cut-off scale of A ~ 0.65 GeV which parameterizes the necessary 
short range dynamics. Most surprisingly, the prediction for the asymmetry energy, A(kf ) = 
33.8 MeV, is in very good agreement with its empirical value. Furthermore, as a nontrivial 
fact pure neutron matter is predicted to be unbound and the corresponding equation of state 
agrees roughly with that of sophisticated many-body calculations for low neutron densities p n < 
0.25 fm~ 3 . In a subsequent work ||, the (complex-valued) momentum and density dependent 
single-particle potential U(p,kf) + iW(p, kf) (i.e. the average nuclear mean field) has been 
calculated in the same framework. It was found that chiral In- and 27r-exchange give rise to a 
potential depth for a nucleon at the bottom of the Fermi sea of U(0, kfo) = —53.2 MeV. This 
value is in very good agreement with the depth of the empirical optical model potential and 
the nuclear shell model potential. 

Given the fact that many properties of nuclear matter can be well described by chiral nN- 
dynamics treated (perturbatively) up to three loop order it is natural to consider in a further 
step finite temperatures T in order to check whether the first-order liquid-gas phase transition 
of nuclear matter is correctly reproduced by this particular dynamics. Such an investigation is 
the subject of the present paper. 

For the relatively low temperatures T < 30 MeV which are of interest in this context one 
can safely neglect effects from thermal pions or thermally excited nucleon-antinucleon pairs. 
As a consequence, nucleons can be treated non-relativistically and the new parameter, the 
temperature T, enters only through the nucleons' thermal occupation probabilities given by a 
Fermi-Dirac distribution. 

In ref . J7J the diagrammatic calculation of the energy density at T = (as a function of the 
particle density p) has been organized in the number of so-called medium insertions. The latter 
is a technical notation for the difference between the in-medium and vacuum nucleon propagator 
(see eq.(3) in ref.fTj). A thermodynamically consistent extension of such an ordering scheme to 
finite temperatures is to relate it directly to the free energy density pF(p,T), since its natural 
thermodynamical variables are the particle density p and the temperature T. In that case 
the free energy density of isospin symmetric nuclear matter consists of a sum of convolution 
integrals of the form, 

poo poo poo 

pF(p,T) = 4 / dp l p 1 K. 1 d(pi) + / d^ / dp 2 K,2d(p 1 )d(p 2 ) 
Jo Jo Jo 
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+ / dpi / dp 2 / dp 3 IC 3 d(p 1 )d(p 2 )d(p 3 ) + pA(p,T), (1) 
Jo Jo Jo 

with the corresponding kernels JCj and the anomalous contribution A(p, T) to be specified.^ 

2 Let us briefly motivate our approach. The standard procedure in field theory is to calculate the grand 
canonical partition function Z or equivalently the thermodynamic potential f2(/x, T) = —(T/V) In Z as a function 
of its natural variables, the (non-relativistic) chemical potential p and the temperature T . In the case of nuclear 
matter at T — the functional relationship £l(p, 0) is however double-valued for p < 0. Since perturbation 
theory usually gives single- valued thermodynamic functions it is more appropriate to calculate first the (single- 
valued) energy density pE (kf) of nuclear matter as a function of the particle density p, as done e.g. in refs.J?], |j] . 
Then eq.(l) defines a thermodynamically consistent extension to finite temperatures T, with the correct T = 
limit. This calculational scheme can also be understood such that in the Legendre-transformation from Q(p, T) 
to the free energy density, p F(p, T) = T) — p dQ/dp, the derivative term is taken into account only for the 
non-interacting (free nucleon gas) part and the (static) one-pion exchange contribution, but not for the higher 
order contributions coming from chiral 27r-exchange. 
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Figure 1: The one-pion exchange Fock-diagram, the iterated one-pion exchange Hartree- and 
Fock- diagrams and the anomalous Fock-diagram. The combinatoric factors of these diagrams 
are 1/2, 1/4, 1/4 and 1/2, in the order shown. 
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denotes the density of nucleon states in momentum space. It is the product of the temperature 
dependent Fermi-Dirac distribution and a kinematical prefactor pj/2ir 2 which has been included 
in d(pj) for convenience. M = 939 MeV stands for the (free) nucleon mass. The particle density 
p is calculated from the density of states in momentum space as 



p = 4 / dpipi d(pi) 
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and this relationship determines the dependence of the effective one-body "chemical potential" 
p(p,T) on the thermodynamical variables p and T. The "true" chemical potential is different 
and given by the formula p = F + p OF /dp. Concerning eqs. (1,2,3) we are following here 
(partially) the approximation scheme of ref. |10|. Our approach is (by construction) thermody- 
namically consistent, since all thermodynamic quantities are derived (via standard relations) 
from the free energy density pF(p,T). The infinite series Li u (x) = Ylk=i k~ u x k defines the 
polylogarithmic function of index v for \x\ < 1. 

The one-body kernel /Q in eq.(l) represents the contribution of the non-interacting nucleon 
gas to the free energy density and it reads 
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While the first two terms are standard |L0 
according to the following criteria. First 



, the correction term, — p\/8M 3 , has been constructed 
it ensures the correct T = limit for the energy per 



particle E(k/) up to order kj, in which p(p,0) = k^/2M and p = 2kj/3Tr 2 . Secondly, the so- 
constructed kernel /Ci combined with the non-relativistic Fermi-Dirac distribution (see eq.(2)) 
gives a very accurate approximation of the free energy density of a fully relativistic free nucleon 
gas 0] for the densities and temperatures of interest here. The factor 4 in eqs. (1,3) counts the 
spin-isospin multiplicity of a nucleon. 



The (non-anomalous) contributions to the free energy density pF(p,T) arising from pion- 
exchange interactions are encoded in the kernels K.2,3 in eq.(l). The closed vacuum diagrams 
related to one-pion exchange (Fock-diagram) and iterated one-pion exchange (Hartree- and 



Fock-diagrams) are shown in Fig. [TJ. The l7r-exchange Fock-diagram including the relativistic 
l/M 2 -correction leads to the following contribution to the two-body kernel /C 2 : 
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As in ref.0, we choose the value g& = 1.3 for the nucleon axial vector coupling constant, /„. = 
92.4 MeV denotes the weak pion decay constant and m n = 135 MeV stands for the (neutral) 
pion mass. The iterated l7r-exchange Hartree-graph (second diagram in Fig. |l]) contributes to 
the two-body kernel /C 2 with the term: 
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In this expression we have omitted the contribution of a linear divergence proportional to 
the momentum cut-off A (see ref. J7|). All such powerlike terms in A are collected in eq. (jTH) . 
Similarly, the iterated l7r-exchange Fock-graph (third diagram in Fig. |l]) gives rise to a two-body 
kernel /C 2 of the form: 
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The additional diagrams of irreducible 27r-exchange (not shown here, but see Fig. 4 in ref.0) 
generate a contribution to the two-body kernel /C 2 via the expression 
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obtained from solving the pion-loop integrals. The complete expression for the power diver- 
gences specific to cut-off regularization reads: 
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The term linear in the cut-off A stems from iterated l7r-exchange with a contribution of the 
Hartree- and Fock diagram in the ratio 4:1. The term quadratic in the cut-off, on the other 
hand, originates from irreducible 27r-exchange. Note that the kernel /Q m ec L-(H3) leads to a 



temperature independent contribution to free energy per particle F(p,T) ~ p. Therefore it is 
fully equivalent to a momentum independent NN-contact interaction. 

Next, we come to the three-body kernel JC 3 which incorporates the temperature and density 
dependent Pauli-blocking effects in intermediate NN-states. The iterated l7r-exchange Hartree- 
graph contributes to the three-body kernel /C3 in the form: 
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and from the iterated l7r-exchange Fock-graph one finds, 
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Note that all integrands in representations of JC 2j3 are odd functions of their respective integra- 
tion variable and therefore one could even drop the absolute magnitude on the lower integration 
limits. We also remind that (non-anomalous) terms involving the product of four Fermi-Dirac 
distributions are effectively absent in second order perturbation theory |TT|. Because of the 
antisymmetry of the accompanying energy denominator under the exchange of two pairs of 
momenta these terms integrate to zero. 

Next, we come to the so-called anomalous contribution A(p, T) in eq.(l). This is a special 
feature at finite temperatures [11, 12] with no counterpart in the calculation of the ground state 
energy density pE(kf) at T = 0. From the "anomalous" Fock-diagram in Fig. 1 we derive the 
following contribution to the free energy per particle of isospin symmetric nuclear matter: 
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with the /i-derivative of the thermodynamical potential due to static l7r-exchange, 
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and the second /i-derivative of the free nucleon gas part, 
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The first term in eq.(14) originates from taking into account the (static) l7r-exchange contri- 
bution in the Legendre-transformation from the thermodynamical potential to the free energy 
density and from the perturbative shift of the chemical potential ft — > ft — Q' l7T (p, T)/Qq(p, T) 
(for details on that procedure, see ref.JT^]). We have explicitly checked that the anomalous con- 



tribution A(p, T) vanishes identically at T = for all densities p. This vanishing, A(p, 0) = 0, 
is an automatic consequence of the Kohn-Luttinger-Ward theorem |12| which of course holds in 
our case since the Fermi-surface and the pion-induced interactions are invariant under spatial 
rotations. Furthermore, it is interesting to observe that the temperature and density dependent 
anomalous contribution A(p, T) vanishes identically in the chiral limit = 0. Note also that 
the anomalous contribution A(p, T) is (formally) of the same order in small momenta as the 
contributions to F(p, T) coming from iterated l7r-exchange. It must therefore not be neglected 
in a consistent and complete calculation. 

Via general thermodynamical relations [|J one finally derives from the free energy per par- 
ticle F(p,T) the pressure P(p, T) and the entropy per particle S(p,T) as: 
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We are now in the position to present numerical results for isospin symmetric nuclear matter at 
finite temperatures. We use consistently the same parameters as in our previous work ]/]]. There, 
our only parameter, the cut-off scale A = 646.3 MeV ~ 7f n , has been fine-tuned to the binding 
energy per particle, —E(kjo) = 15.26 MeV. Fig. 2 shows the free energy per particle F(p,T) as 
a function of the nucleon density p for various temperatures T = 0, 10, 18, 25.5, 30 MeV. The 
uppermost line is the attractive branch of the nuclear matter saturation curve at T = 0. The 
singular behavior of the free energy per particle F(p, T) for p — > at T > is a generic feature 
(see e.g. the figures and tables corresponding to the results of the Urbana group in ref.p). The 
internal energy per particle F(p,T) + T S(p,T), on the other hand, approaches approximately 
the value 3T/2 for p — > 0. 

Fig. 3 shows the calculated pressure isotherms P(p, T) of isospin symmetric nuclear matter. 
As it should be these curves display a first-order liquid-gas phase transition similar to that of 
the van-der-Waals gas. The coexistence region between the liquid and the gas phase (which has 
to be determined by the Maxwell construction |4[]) terminates at the critical temperature T c . 
From there on the pressure isotherms P(p, T) grow monotonically with the particle density p. 
In the present calculation we find a critical temperature of T c = 25.5 MeV and a critical density 
of p c = 0.09 fm -3 ~ 0.5po (with p — 0.178 fm -3 , the predicted nuclear matter saturation 
density). Together with the critical pressure P(p c ,T c ) = 0.69 MeV fm~ 3 this prediction of the 
critical point deviates considerably from most other calculations. The too high value of the 
critical temperature T c = 25.5 MeV finds its explanation in the strong momentum dependence 
of the underlying single-particle potential U (p, k/o) near the Fermi-surface p = kfo (see Fig. 3 in 
ref. |Q). The nominal value of the effective nucleon mass at the Fermi-surface is M*(kfo) — 3M. 
The latter quantity determines via the reciprocally related density of states at the Fermi- 
surface completely the low-temperature behavior of a Fermi-liquid. Since in our calculation 
the density of thermally excitable quasi-particles is too low nuclear matter has to be heated up 
to higher temperatures until it evaporates completely. More elaborate calculations of nuclear 
matter in effective (chiral) field theory are necessary in order to cure this particular problem 
of the large effective nucleon mass M*(kj ). Note that despite its large isospin factor 18 the 
anomalous contribution A(p, T) does practically not influence the behavior of nuclear matter 
at low temperatures T < 30 MeV. 

As a side remark, we consider the chiral limit m n = 0. With a reduced cut-off scale of A = 
555.8 MeV (and fixed gA, fn, M) the same maximum binding energy per particle (15.26 MeV) is 
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Figure 2: The free energy per particle of isospin symmetric nuclear matter F(p,T) versus the 
nucleon density p. Each curve is labeled with its corresponding constant temperature T. 

obtained at a saturation density of p — 0.145 fm~ 3 . Interestingly, the critical temperature T c ~ 
27MeV remains nearly unchanged when taking the chiral limit. This confirms the expectation 
that the critical temperature T c is primarily determined by the binding energy per particle 
—E(kfo) and the effective nucleon mass at the Fermi-surface M*(kfo). 

The present calculation is readily adapted to the case of pure neutron matter at finite 
temperatures T. Only the isospin factors of the pion-exchange diagrams (see Fig. 1) change 
when switching from isospin symmetric nuclear matter to pure neutron matter. The free energy 
per neutron F n (p n ,T) takes a form analogous to eq.(l), namely: 



The neutron density p n = p/2 is reduced by a factor 2 in comparison to p calculated via eq.(3) 
and the same applies to the one-body term proportional to /Q in eq.(18). The contributions 
from l7r-exchange and iterated l7r-exchange to the neutron kernels K, n j are determined by 
certain relative isospin factors as, 
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while irreducible 27r-exchange leads to the expression 
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Figure 3: Pressure isotherms P(p, T) of isospin symmetric nuclear matter. The coexistence 
region of the liquid and gas phase ends at the critical point p c ~ 0.09 fm -3 and T c ~ 25.5 MeV. 
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The expression for the power divergences specific to cut-off regularization changes accordingly 
in the case of pure neutron matter, 
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Furthermore, for pure neutron matter the anomalous contribution gets reduced by a relative 
isospin factor 1/18, or written in an explicit formula: A n (p n) T) = A(2p n ,T)/9 with A(p, T) 
given in eq.(14). With this small weight factor the anomalous 27r-exchange Fock diagram 
becomes irrelevant for pure neutron matter. 

In Fig. 4, we show by the full lines the calculated free energy per neutron F n (p n , T) for tem- 
peratures T = 10, 20 MeV. The dashed lines in Fig. 4 correspond to the many-body calculation 
of the Urbana group 0. These curves should be considered as a representative of the host 
of existing neutron matter calculations which scatter around them. In order to demonstrate 
the effects of the nn-interaction, we show by the dotted lines in Fig. 4 the result of a free neu- 
tron gas (as given by the one-body kernel /Q in eq.(18)). One observes a fair agreement of 
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Figure 4: Free energy per particle of pure neutron matter F n (p n ,T). The full lines give the 
result of chiral one- and two-pion exchange. The dashed lines correspond to the many-body 
calculation of ref . || and the dotted lines show the result for a non-interacting neutron gas. 



our calculation with the results of ref.0 for low neutron densities, p n < 0.2 fm -3 . One may 
also conclude from Fig. 4 that the parameterfree interaction effects generated by chiral In- and 
27r-exchange are fairly realistic in this low density region. 

In summary, we have extended our recent three-loop calculation of nuclear matter in chi- 
ral perturbation theory to finite temperatures T. The contributions from one- and two-pion 
exchange diagrams are included in the density and temperature dependent free energy per 
particle F(p,T). This guarantees thermodynamical consistency and the correct T = limit 
to the energy per particle E(kf). We reproduce the familiar liquid-gas phase transition of 
isospin symmetric nuclear matter. The predicted critical point, T c = 25.5 MeV, p c = 0.09 fm~ 3 
and P(p c ,T c ) = 0.69 MeV fm -3 lies however too high in temperature. The obvious reason for 
that is the strong momentum dependence of the single-particle potential U(p, kfo) || near the 
Fermi-surface p = k/o implying a huge effective nucleon mass M*(kfo) ~ 3M. More elaborate 
calculations of nuclear matter in effective (chiral) field theory are necessary in order to cure this 
problem. Furthermore, we have considered pure neutron matter at finite temperatures T in the 
same framework. We have found fair agreement with sophisticated many-body calculations for 
the low neutron densities, p n < 0.2 fm -3 , relevant for conventional nuclear physics. 
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